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Abstract 

The mid-rapidity charged particle multiplicities in pp and AA collisions at LHC energies are described 
in the framework of a generalized eikonal model with shadowing corrections incorporated in AA. We show 
that the pp data require a Pomeron intercept close to 1.2, higher than the conventional one close to 1.1. An 
s°'^^ energy dependence is obtained in the LHC range and beyond. The size and centrality dependence of 
the AA multiplicity at y/s — 2.76 TeV is reproduced and its energy dependence is predicted. 



1 Introduction 

Recent data on mid-rapidity charged particle multiplicities in proton-proton collisions by the ALICE [1] [2] , 
CMS |3] and ATLAS |1] collaborations have extended the energy range of this observable up to y/s — 7 TeV. 
They show an s-dependence in s°-^^. The same dependence is observed at lower energies. Such an s-dependence, 
observed at high energies, has important consequences for the Pomeron intercept. This is based on a very general 
property of multiple scattering models, known as the AGK cancellation |5]. As a result of this cancellation, 
absorptive corrections (i.e. multiple-scattering effects) vanish identically in the single particle inclusive cross- 
section da/drj. This theorem is valid for a large class of multiple-scattering models, including eikonal and 
Glauber type models. It is the equivalent in soft processes of the factorization theorem in pcrturbative QCD. 

Consider the charged multiplicity dN^^/drj — {da^^ / drj) / apj^ . The data show an s°'^^ behaviour for 
dNPP/dr]. Assuming that the total non-diffractive pp cross-section a^J^ behaves as s"" with a^. ^ 0.07 — 0.09, 
as observed at energies below LHC, we conclude that da^^ /d-q behaves as with A ^ 0.18 — 0.20. Due to 
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the AGK cancellation this is the behaviour of single scattering (Born term) and thus the Pomeron intercept is 



String models like DPM ^ or QGSM [5] do, of course, satisfy the AGK cancellation. This cancellation is fully 
effective at mid-rapidities and high energies. Energy conservation introduces violations in the fragmentation 
regions which, at low energies, propagate to mid rapidities. It turns out that some amount of violation is still 
present in mid-rapidity AuAu collisions at RHIC which go away at mid-rapidities in the LHC energy range. 

Consider next processes involving nuclei. Here the AGK cancellation implies that the charged multiplicities 
dN^^/drj scale with the number of binary collisions, i.e. the same scaling behaviour valid in hard processes 
due to the factorization theorem. It is well-known that such a scaling is not supported by the data. Actually, 
the recent ALICE data [11], lately confirmed by CMS il2] and ATLAS [13^ collaborations, show that the 
mid-rapidity charged multiplicity per participant pair — dN^^ / drj/{npart/'^) in central PbPb collisions 
at ^/s — 2.76 TeV is about two times larger than in pp collisions at the same energy. However, the AGK 
cancellation implies that this ratio should be equal to the ratio ricoii / {ripart / 2) which is about five at RHIC and 
increases with energy. The factor of two in Ra relative to pp, conveys the idea of a greater efficiency of AA 
collisions for particle production in each binary collision. However, the surprising fact is the strong suppression 
of the AA multiplicity as compared to the predictions of the eikonal- Glauber model. The solution of this puzzle 
is the presence in AA collisions of the so-called shadowing corrections, not included in that model. These 
corrections are small in pp but are quite large in AA collisions where they are enhanced by A^^^ factors. In 
view of that, it is clear that an analysis of the results mentioned above puts strong constraints on the size and 
energy dependence of shadowing. 

The plan of the paper is as follows. In Section 2 we describe the physical bases of both AGK cancellation 
and shadowing. In Sections 3 and 4 we study the mid-rapidity pp and AA charged multiplicities, respectively. 
Conclusions are given in Section 5. 

2 AGK cancellation and shadowing 

Let us recall the physical bases of the AGK cancellation ^ and shadowing [101 E] • 

A simple way to illustrate the AGK cancellation is to consider the cross-section af/^{h) for v inelastic collisions 
in pA scattering in the probabilistic Glauber model 



This formula is self-explanatory. Ta(6) is the nuclear profile function normalized to unity. The first factor in 
eq. ([l]) is the number of ways to choose v interacting nucleons out of A. The second one in the probability 
that V nucleons interact at fixed impact parameter h. The third factor is the probability for no interaction of 
the remaining A — v nucleons. Upon summation of ([T]) from v = 1 to A and integration in 6 one obtains the 
non-diffractive pA cross-section. It can be seen numerically that it behaves like A" with a ^ 2/3. On the 

tA similar pomeron intercept, ap(0) = 0.2 fg], has been prevoiusly used by ALICE collaboration [J] to fit the central rapidity 
p/p ratio as a function of the center-of-mass energy. 



ap(0) = 1 + A|t| 




(1) 



3 



contrary for the multiplicity we have 

^ — y ,yaP/ ^ = V ^ (2) 
di] GpA ^ drj drj 

with V = Affppl GpA- The factor v in eq. ([2]) is due to the fact that, in the presence of v non-diffractive inelastic 
collisions, the trigger particle can be produced in any of them. Eq. ([l]) is equivalent to da^^/di] — A daP^/drj. 
Thus dcjP^ / drj scales with A and dNP^/drj scales with the number of binary collisions. In the former observable 
all multiple scattering contributions cancel identically and one is left with the scaling in A^ of the Born term. 
The same cancellation takes place in AA collisions in the Glauber-Gribov model (see Section 4) and in pp in 
the eikonal model (see Section 3). In the framework of Reggeon Field Theory [IS] the AGK cancellation is 
illustrated in Fig. [l] for a double scattering. 




a) b) 

Figure 1: A double scattering (two-Pomeron exchange graph) in pA collisions with a single inelastic collision 
(one cut Pomeron), a), and two inelastic collisions (two cut Pomerons), b). 

The first graph, with one cut Pomeron (i^ ~ 1), corresponding to a single inelastic collision, has a weight — 4C 

A 

while the second one, corresponding to two inelastic collisions, has a weight +2C. In the summation ^ a^,-^ 

v=l 

A 

there is a net contribution — 2C (absorptive correction) which reduces the A dependence, while in ^ (jP/^ 

v=\ 

the net contribution is zero. This result is true for any multiple scattering graph and for a very general class 
of blobs in the upper part of the diagram. This is the origin of the AGK cancellation [5] in Reggeon Field 
Theory [H]. 

Let us turn next to shadowing effects. They are related to graphs associated with large-mass diffraction, 
commonly known as triple Pomeron graphs, as the ones represented in Fig. [2] In Fig. [2^ the trigger particle is 
produced at a laboratory rapidity smaller than the one of the junction of the three Pomerons. In this case we 
are in the same situation as in the graphs of Fig. [TJand the AGK cancellation takes place. The contribution of 
this type of graphs to da^^ /dr] cancels identically and we recover an A^ behaviour. On the contrary, in the case 



a) b) 

Figure 2: A triple Pomeron graph for daP^ /drj with the laboratory rapidity of the produced particle smaller 
(a) or larger (b) than the one of the triple Pomeron vertex. 

of Fig.[2]D, the trigger rapidity is larger than the one of the triple Pomeron vertex and the factor v is not present. 
In this case absorptive corrections are at play which reduce both the size and the A dependence of a^^ /drj. 
The contribution of the triple Pomeron graphs in Fig. [2] can be computed using parameters determined from 
large-mass single diffraction data. Following refs. [TUl [S] we use Schwimmer unitarization scheme [T^]. The 
suppression from shadowing in AA collisions for a particle produced at mid-rapidity {y* = 0) is then obtained 
replacing the nuclear profile function Taa(^) = / <P sTA{s)TA{b — s) by 

^ '"'~J l + AF{y*=0)TA{s) l + AF{y*=0)TA{b~s) 

where 

F{y* = 0) = C [exp{Ay^ax) - exp(Ay™„] /A (4) 

is the triple Pomeron graph contribution, with the rapidity of the triple Pomeron vertex integrated up to y* = 
(where the trigger particle is produced), i.e. up to ymax — U* + in{^/s/mT)- We take ymin = in{RA'mN / VS)- 
Smaller values of y, corresponding to larger values of the mass of the diffractively produced system, are cut-off 
by the nuclear form-factor {tmin effect). Here tot is the transverse mass of the particle[^| rnjv the nuclear mass 
and Ra = 0.82 A^^^ + 0.58 fm - the gaussian nuclear radius. TA{b) is the nuclear profile function for which we 
use a Woods-Saxon parameterization [17] . 

Let us discuss now the values of A and C. In previous works [T3] we have used A = 0.13 and C = 0.04 fm^ 
(C/A = 0.31 fm^) corresponding to the Pomeron intercept ap(0) = 1.13. As discussed in the Introduction we 
show in Section 3 that the LHC data require a higher Pomeron intercept ap(0) = 1.19. Consequently we take 
A — 0.19. Note that we have already used a similar value of A in 18j to describe 7*p data, A = 0.2. With this 
value of A introduced in the triple Pomeron a good description of the M^-dependence of high-mass diffraction 



^Wc take the transverse mass of the pion, my = 0.38 + 0.0233 * hi{y/s/5^) GeV. 
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production at = was obtainecj^ This increase in A leads to a decrease in the residue of the Pomeron -to 
which the constant C is proportional. This induces a reduction of its value to C = 0.03 fm^. With these values 
of C and A the values of the shadowing suppression factor differ from the old ones by less than 10% in the 
energy range between y/s — 200 GeV and 5.5 TeV. 

We see from eq. (jsj that shadowing corrections are enhanced in nuclear collisions by factors A^/^^, and for 
this reason, they turn out to be quite important in this case. They are small in pp and will be neglected in the 
calculations developed in Section 3. 



3 Multiplicities in pp collisions 



Multiple scattering in hadron-hadron collisions is usually described in a generalized eikonal model in which 
inelastic intermediate states are included in the vertex functions of the multiple scattering graphs. This cor- 
responds to the inclusion of diffractively produced intermediate states. In Reggeon field theory the large-mass 
diffractive states correspond to triple Pomeron (see Fig. 2) and loop graphs. A description of this approach, 
with exponential residues in t -corresponding to gaussians in impact parameter- can be found for instance in 
[T5] and in Appendix B of ref . [5] . A simpler albeit cruder approach is obtained by neglecting the contributions 
involving triple Pomeron couplings and including low-mass diffractive contributions as intermediate states in 
the eikonal model [H1[5T]. Using exponential residues in t, and a Regge behaviour for single scattering (Born 
term) in (s/sq)"^'*''"^ with the Pomeron intercept ap{t) = 1 + A + a'pt, all the loop integrals can be performed 
analytically. The cross-sections (7^ for k non-diffractive inelastic collisions -the equivalent of (jP-^ in eq. ([l])- 
are then given by [SI |5T] 



Op 

kZ 



l-exp(-Z)^ 



(5) 



Here ^ = £n{s/so) with Sq = 1 GeV^, crp = 87r7p exp(A^), and Z 
From ([5]) we obtain the non-diffractive inelastic cross-section 

K>1 

The mid-rapidities dN^P /drj is proportional to k = ^ k ctk/ ^ CTk, as we will explain below. 



(6) 



An important parameter in ^ is the Pomeron intercept. We take A = 0.19. Such a large value has been 
motivated in the Introduction. We also take a'p — 0.25 GeV^^ and — 3.3 GeV^^ OET]. These parameters 
control the t-dependence of the elastic peak and its energy dependence (shrinking). With these values for 
and a'p one obtains [S] a value of the elastic pic slope B ~ 20 GeV~^ at y/s = 7 TeV, in very good agreement 
with the recent TOTEM measurement [22]. Note that this value of B is larger than the one in the Born term 
due to the effect of the unitarity corrections. The parameter jp = 0.85 GeV^^ controls the size of the interaction 
cross section and has been determined from the absolute normalization of cr^^ , which leads to the values given 
in Table [T] Note that the energy dependence of a^^^ is essentially determined by the Pomeron intercept A, 

§Notc that an (l/M^)^-^ dependence for large-mass diffraction has been obtained |19l from a fit of the data. However, it has 
been shown [20) that the power of is larger in the non-absorbed (triple Pomeron) graph than in the absorbed one, i.e. the 

one measured in experiment. 
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and it is about s^ *^^, in agreement with experimental data. Finally, the quantity Ce takes iirto account the 
modification of the eikonal due to inelastic diffractive states. In refs. 9 and 2 1 a value Ce = 1.5 is used 
corresponding to a 50% contribution of low-mass diffractive states relative to the elastic ones. Note that here 
we have also included in Ce high-mass diffraction and taken Ce = 1.8. Actually, the percentage contribution 
of large-mass diffractive states, given by the triple Pomeron graph, increases with energy. However, this energy 
dependence turns out to be small in the LHC energy rang^ Note also that including the large-mass diffraction 
contribution in C^; is only possible if the ^-slopes of the two-Pomeron-exchange graph and the triple Pomeron 
graph are equaji[[ Diffractive cuts corresponding to fc = have not been included. A treatment of diffractiction 
would require a more refined approach, proceeding for instance along the lines of refs. [181 123| . 

In string models the proportionality between the pp multiplicities and k mentioned above is expressed as [8] 
dNPP/dr] — 2KdNs /di], where the last factor is the multiplicity in a single string. This relation is valid at high 
energies -when all strings have equal mid-rapidity multiplicity. dNs/dr]{y* = 0) is a constant, independent of 
s. In the calculations we have used 2dNs/dri = obtained from the experimental value dN^P/drj = 4.7 at 

= 2.76 TeV [1]. The results for and dNPP/d'q{y* = 0) as a function of s are given in Table[l]and Fig.js] 




2 3 4 5 

10 10 10 10 

s^'^(GeV) 

Figure 3: The values of dNP^/dr] obtained from the formulae in Section 3 are plotted as a function of y/s in 
the range ^/s = 100 GeV to 50 TeV and compared to non-single-diffractive (NSD) data [H [3l [24l [25l [26 l [27 ] . 
The line s"'^^ (dashed line) representing the available experimental data is also plotted. 

We see that in the LHC energy region and beyond the pp multiplicity follows closely the s^'^^ behaviour 
observed in the data |1,. Actually, this energy behaviour is predicted to hold in a much larger energy range 
-at least up to 50 TeV. As a consequence cTpJ^ follows closely an s° °^ energy dependence. Note that at lower 

lllntegrating the dependence of the triple Pomeron [M^]~''-~^ from = 5 GeV^ to s/20, with A = 0.19 we get a 6% 
increase between -y^ = 2.76 and 14 TeV. 

lit should also be mentioned that with A ^ the signature factor induces a well defined real part in the scattering amplitude. 
However, only its imaginary part contributes to the non-diffractive observables considered in this paper. 

**This is to be compared with a value of 2dNg /dy{y* = 0) about 2 used in DPM, corresponding to a plateau height of a single 
string close to one [10) . Note that the ratio dN^^ /dy over dN^P /drj is close to 1.3 at y* = 0. 
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Table 1: The values of the non-diffractive pp cross-section and charged particle pseudo-rapidity densities in the 
central rapidity region for the energy range between 200 and 50000 GeV. 



y/s (GeV) dNPP/dr]{y* = 0) aE," (mb) 



200 


2.99 


31.22 


540 


3.50 


38.97 


900 


3.82 


43.33 


1800 


4.34 


49.64 


2760 


4.71 


53.77 


5500 


5.42 


60.78 


7000 


5.70 


63.33 


14000 


6.61 


70.99 


50000 


8.80 


86.19 



energies the calculated multiplicity is larger than the measured one. This is consistent with energy-momentum 
conservation effects which reduce the former one. At ^/s ~ 200 GeV the difference between calculated and 
observed multiplicity is about 15%. 

Note that our calculated values of dN^P/dy refer to non-diffractive (ND) multiplicities. However, the avail- 
able experimental data on multiplicities refer to non-single-diffractive (NSD) interactions, or to inelastic (INEL) 
ones. Taking into account the smallness of the double-diffractive contribution (DD), [212111511121, the difference 
bewteen non-diffractive (ND) and non-single-diffractive (NSD) multiplicities will be small. Therefore, we have 
compared our theoretical non-diffractive value (ND) to experimental non-single-diffractive (NSD) data. 



4 Multiplicities in PbPb collisions 

The AGK cancellation implies that at mid-rapidities da^^/drj ~ A^da^P/drj, which implies 

dN^^ „ cr^^ dNPP dNPP 

^ ^2 ^ ^ ^^^^^ 

drj aAA drj drj 

where ricoU = '^pjF I '^aa is the average number of binary nucleon-nuclcon collisions. As a function of impact 
parameter we have 

ND 

n,M{b) ^ ^^TAA{b) (8) 
caaKo) 

with aAA{b) = 1 - exp[-cr^^A2 TAA{h)] and TAA{b) = J d'^s Ta{s) TAib - s) (see Section 2). It follows from 
eqs. ([T]) and (|8| that the multiplicity per participant pair is given by 

dr] npartip) drj 

As discussed in the Introduction, the first factor in eq. ^ is about 5 at RHIC for central AuAu collisions 
and increases with s. This is much larger than the factor 2 observed experimentally. However, as discussed in 
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Section 2, the AGK cancellation is only valid in the absence of shadowing. In its presence, eq. Q has to be 
modified. One has 

^ib)^A'^^S^\b) (10) 
dv (^AA(b) dr] 

where the shadowing correction S'''''(&) is given by eqs. ^ and In this way, the AA multiplicity is strongly 
reduced [TOllM]. This reduction is given by the suppression factor S'^^{h) /TAA{b)- Its numerical values for 
PbPb collisions at ^/s = 2.76 and 5.5 TeV are given in Table [2] as a function of h. 




Figure 4: The values of dN^^/dr] at y/s — 2.76 and 5.5 TeV are plotted for different centrality bins. The 
dashed lines are obtained from eq. ( 10 ) with shadowing corrections given by eqs. (|3| andQ and using the 



non-diffractive values dN^^/dr] — 4.7 at -y/s = 2.76 GeV and 5.4 at y/s — 5.5 TeV. The full lines are obtained 
using instead the inelastic values dN^P /drj — 3.9 at = 2.76 GeV and 4.5 at = 5.5 TeV (see main text). 
Data are from [H] (open symbols) and 13J (full symbols). The values for inelastic and non-single-diffractive 
dNPP/drj at = 2.76 TeV are also shown 0. 

Using the numerical values in Tables [l] and [2] we can compute the mid-rapidity AA multiplicities for any 
value of the impact parameter. The results at y^s = 2.76 and 5.5 TeV are shown in Fig. 4 (dashed lines). At 
y/s ~ 2.76 TeV the centrality dependence of the data [5] is well reproduced. Practically the same centrality 
dependence is predicted at ^/s = 5.5 TeV. The same feature has been observed experimentally between ^/s = 
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Table 2: The values of the shadowing suppression factor S'^{h)/TAA{h) for PbPb cohisions at = 2.76 and 
5.5 TeV, as a fimction of the impact parameter 6, computed from eqs. (Is]) and Q. 



b = 2.76 TcV ^fs ^ 5.5 TeV 



0.2828 0.2390 

2 0.2839 0.2396 

4 0.2910 0.2467 

6 0.3091 0.2637 

8 0.3432 0.2959 

10 0.4091 0.3599 

13 0.6167 0.5714 



200 GeV and 2.76 TeV [2 . Its increase between ^/s = 2.76 and 5.5 TeV is slightly less than the corresponding 
increase in the pp multiplicity. An increase larger than the one in pp has been observed between ^/s = 200 GcV 
and 2.76 TeV. 

The absolute value of the multiplicity at ^/s — 2.76 TeV is about 25% higher than the experimental one. 
Note, however, that there is an uncertainty in this absolute value due mainly to the uncertainty in cr^J^. 
Moreover, there is a puzzling feature in the data |11| . namely the ratio between the two extreme centrality 
bins multiplicities is smaller than the one between the most central bin and pp. Therefore the multiplicity in 
peripheral bins gets smaller than the pp one. Actually the data extrapolate nicely to the inelastic pp multiplicity 
(3.9 at y/s = 2.76 TeV) rather than to the non-diffractive one (4.7 at y/s = 2.76 TeV). Using the former value 
in our calculations, instead of the non-diffractive one, we get absolute values in agreement with experiment. 

5 Conclusions 

We have studied the mid-rapidity charged particle multiplicities in pp and AA collisions in the LHC energy 
range. We have shown that the observed energy dependence of the pp multiplicity requires a Pomeron intercept 
close to 1.2. One obtains in this way an energy dependence remarkably close to s*^'^^, in good agreement with 
available data. This behaviour is predicted to hold up to at least 50 TeV. At lower energies the calculated pp 
multiplicity is larger than the experimental one. The difference is about 15 % at -y/s = 200 GeV. This behaviour 
is consistent with energy conservation effects which reduce the multiplicity and are present at low energies. 
These effects go away when we reach the LHC range. 

In AA collisions the centrality dependence of the multiplicity is predicted to be practically identical at 
•y/fi = 2.76 and 5.5 TeV, in agreement with lower energy extrapolations. Its size is well reproduced. The 
increase in multiplicity with energy between 2.76 and 5.5 TeV is predicted to be slightly less than the one in 
pp, while a larger increase is expected from lower energy results which show an energy dependence in s^'^^ for 
AA and " for [H]. 

The obtained Pomeron intercept is larger than the conventional one used in most models obtained from lower 
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energy data. In this respect one can consider two different scenarios. A first possibility 6 is that a Pomeron 
intercept close to 1.2 can describe the data at all energies when all multiple scattering contributions are taken 
into account -together with low energy corrections such as energy conservation effects in the multiplicities. A 
second scenario is the so-called Pomeron flavoring |30| . In this scenario the Pomeron intercept increases due 
to the opening up of effective thresholds for heavy flavor production as the energy increases. Actually the 
incorporation of semi-hard events (mini-jets) within DPM |31|, I32j can be considered a flavoring phenomenon 

m- 

Calculations of the pp multiplicity at LHC in the framework of the QGSM [5] can be found in [34 . In these 
calculations the Pomeron intercept is taken to be 1.12. As a consequence, the energy dependence of the pp 
multiplicity is substantially smaller than ours. Monte Carlo models such as PYTHIA [35] and PHOJET [36] also 
have a too small energy dependence in the LHC range [3]. All these results lend support to the necessity of a 
higher Pomeron intercept. Obviously our value of an intercept close to 1.2 has to be validated by a measurement 
of the energy dependence up to 14 TeV. 
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